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ABSTRACT 

We present deep infrared (Ks band) imaging polarimetry and radio (1.4 and 4.8 GHz) po- 
larimetry of the enigmatic transient Swift J164449. 3+573451. This source appears to be a 
short lived jet phenomenon in a galaxy at redshift z = 0.354, activated by a sudden mass 
accretion onto the central massive black hole, possibly caused by the tidal disruption of a 
star. We aim to find evidence for this scenario through linear polarimetry, as linear polari- 
sation is a sensitive probe of jet physics, source geometry and the various mechanisms giv- 
ing rise to the observed radiation. We find a formal Ks band polarisation measurement of 
Pii n = 7.4 ± 3.5% (including systematic errors). Our radio observations show continuing 
brightening of the source, which allows sensitive searches for linear polarisation as a func- 
tion of time. We find no evidence of linear polarisation at radio wavelengths of 1.4 GHz 
and 4.8 GHz at any epoch, with the most sensitive 3cr limits as deep as 2.1%. These up- 
per limits are in agreement with expectations from scenarios in which the radio emission is 
produced by the interaction of a relativistic jet with a dense circumsource medium. We fur- 
ther demonstrate how the polarisation properties can be used to derive properties of the jet in 
Swift J164449. 3+573451, exploiting the similarities between this source and the afterglows 
of gamma-ray bursts. 

Key words: gamma-ray burst: individual: GRB1 10328, techniques: polarimetric, galaxies: 
jets 



1 INTRODUCTION 

Recently, the Swift satellite triggered on a peculiar new gamma-ray 
source. While initially thought to be an unusual gamma-ray burst 
(GRB 110328A; Cummings et al. 2011), prior detections by the 
Burst Alert Telescope (BAT) in addition to X-ray behaviour highly 
atypical for GRBs (see e.g. Burrows et al. 201 1; Levan et al. 201 1) 
showed this source was likely of a different nature. Localisation of 
a variable infrared and radio counterpart (Levan et al. 201 1; Zaud- 
erer et al. 201 1) spatially coincident with the nucleus of a galaxy at 
redshift z = 0.3543 (Levan et al. 2011), further motivated a model 
in which the transient represents a sudden activation of a blazar-like 
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phenomenon, most likely fed by the tidal disruption of a (main se- 
quence) star l Bloom et alj|201lt ICannizzo. Troia & Lodatdl201ll : 



iKrolik & Pirarj|201 ll) . though it is difficult to rule out a completely 
novel phenomenon more closely related to standard gamma-ray 
bursts JOuataert & Kasenll201ll ; lwooslev & Hegej|201ll) . 

Whereas the gamma-ray event itself, as well as the early X- 
ray behaviour of the transient, is not in accord with predictions of 
tidal disruption events, the general behavio ur at later times se ems 
to match the expectations for these events l lBloom et alj|201 fh . In 
this model, the sudden accretion onto the central massive black 
hole caused by the tidally disrupted star gives rise to a relativistic, 
beamed outflow whos e light dominates o ver the light produced by 
the stellar disruption (Bloom et al. 201 1). The spectra l energy dis- 
tribut i on of the transient shows two prominen t bumps dBloom et al.l 
1201 ll ; IZaudereJ l201ll ; iBurrows et all l201lh . which may be ex- 
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plained by the presence of synchrotron and inverse Compton emis- 
sion p rocesses dBloom et alj|201 ll ; iBurrows et alj|201 ll ; lAliu et al 



2011 ). Constraints on the source size from radio spectra dZauderei 



20111) imply a relativistic jet as the source of the observed syn- 



chrotron emissio n, in agreement with fits to the source spectral en- 
ergy distribution (Bloo m et al.l201 ll ; IZauderej201 lklBurrows et al.l 
1201 lb . This emission is likely generated in the shock interaction be- 
tween the jet and the circumsource medium, in contrast to blazars 
in which the synchrotron emission is produced within the jet it- 
self. Modelling of the radio emission indicates a start time of the 
initial event some 4 days prior to the time BAT triggered, in agree- 
ment with analysis of BAT data take n before the trigger occurred 
l lBurrows et al.l(201ll ; lza"uderej|201lb . 

The linear polarisation of jet sources gives crucial insight into 
the acceleration mechanism of the jet, the configuration of the mag- 
netic field responsible for the synchrotron emission, and the overall 
geometry of the jet with respect to the line of sight. In the case 
of this source it may enable us to establish the position of Swift 
J164449. 3+573451 within the family of accreting jet sources, in 
particular the relation with blazars and X-ray binaries. In this paper 
we study the linear polarisation of Swift 1164449.3+573451 using 
deep Ks band imaging polarimetry and radio polarimetry, frequen- 
cies at which dust obscuration has the smallest influence and at 
which other transient jet sources are thoroughly studied. 

This paper is organised as follows: in Section 2 we detail our 
observing strategy, data reduction and analysis; in Section 3 we 
compare the polarimetry with models for Swift J 164449.3+573451 
and other jet sources; and in Section 4 we present our conclusions. 



2 OBSERVATIONS AND DATA ANALYSIS 
2.1 WHT LIRIS Ks band polarimetry 

Levan et al. (2011) show that in the optical range of the spectrum 
the transient is highly extinguished by dust in the host, and large 
scale variability is therefore only evident in the near-infrared wave- 
lengths. This makes optical polarimetry of the transient unfeasible. 
Infrared polarimetry at the longest wavelengths (i.e. K band) is pre- 
ferred, despite the greater difficulty in obtaining polarimetry with 
sufficiently small errors of faint sources in this wavelength range. 

We acquired deep Ks band polarimetry of the source using the 
4.2m William Herschel Telescope at La Palma, Spain, starting on 
02:13 UT on April 14th 2011. We used the LIRIS (Long-slit Inter- 
mediate Resolution Infrared Spectrograph; Manchado et al. 2004) 
instrument in imaging polarimetry mode, which utilises a double 
Wollaston configuration (see Oliva 1997) to obtain simultaneous 
measurements of the polarised flux at angles 0, 90, 45 and 135 de- 
grees - note that LIRIS does not possess a half-wavelength plate at 
the time the data were acquired. An aperture mask is used to avoid 
overlap of the different angle images, resulting in a field of view of 
1x4 arcmin. 

We express the polarimetry information in terms of the Stokes 
vector? = (Q, U, V,I) (see e.g. Chandrasekhar 1950), where the 
components of this vector have the following meaning: Q and U 
contain the behaviour of the linear polarisation; V the circular po- 
larisation and / is the total source intensity. Generally we will use 
the normalised Stokes parameters Q/I and U/I in this paper. The 
LIRIS setup means that these can be found from 

aw / C0-/90) , .... (/45-J135) 
1 = Q/I = — ; — 7 and u = U/I ~ 




Figure 1. A combined image of a subset of the data of rotl (-90 degrees 
position angle), with the transient indicated with tick marks. 



these Stokes parameters can be determined from one single expo- 
sure, under equal conditions, which is particularly advantageous 
at infrared wavelengths where sky conditions can change rapidly. 
Theoretical models are often expressed in terms of the linear po- 
larisation degree (or fraction) P\ m and polarisation angle 9. These 
quantities can be found from the Stokes parameters through: 

1 
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(Us + /135) 



This shows clearly the advantages of a double Wollaston design: 



Py m = ^jq 2 + u 2 and 8 = -arctan(;(/g). 

Note that the conversion from Stokes parameters to P lm brings with 
it complications, discussed further below, so wherever possible we 
will work in Stokes parameter space. 

We obtained a total of 7340 sec exposure time, using 20 sec 
integrations and a 5 point dither pattern under fair conditions (av- 
erage Ks band seeing 0.8 arcsec). The dithering was primarily in 
the mask X direction so the transient never disappeared behind the 
aperture mask. The field orientation, -90 degrees, was chosen such 
that several bright foreground stars fall within the 1x4 arcmin aper- 
ture mask (Fig.[T](, so that they can serve as secondary calibrators. 
Approximately half the observation, total 3400 sec exposure time, 
of which three exposures could not be used due to trailing of the 
stars, were taken with -90 degrees rotation. The remainder of the 
observations were done with field rotated a further 1 80 degrees to 
decrease the effects of flatfielding errors and imperfect behaviour of 
the Wollaston prisms: a 180 degree instrument rotation corresponds 
to a full rotation in q, u space. Hereafter we will call the datasets 
with the two different rotation angles rotl and rot2, respectively. 

To calibrate the resulting polarimetry we used observations of 
zero polarisation standard star BD +33 2642. Observations of this 
standard consisted of 5 exposures of each 5 sec, starting at 01:21 
UT on April 15th 201 1, and were taken without instrument rotation 
(i.e North up, East left). 

We clipped from the raw frames the images corresponding to 
the four different angles, and treat these separately but in an iden- 
tical way. We flat field corrected all data using tasks within IRAF 
using twilight sky flats taken the same night. We then performed 
sky subtraction by constructing sky frames per 5 dither cycle. The 
resulting sky subtracted frames were aligned and combined using 
our own software, separately for the rotl and rot2 datasets. In the 
same manner we processed the data taken of the zero polarisation 
standard star. 

We used our own software in combination with IRAF routines 
to perform aperture photometry on the images of the transient to 
find the source fluxes / , 7 90 , /135 and / 45 . We used a seeing matched 
aperture of 5 pixels radius (1.5 times the on-frame full width at 
half maximum (FWHM) of the point spread function). A Stokes / 
image was created as well; photometry performed on this image is 
shown in the lightcurves in Fig. [3] It is clear that polarimetry of an 
object this faint is highly challenging. 

We compute the Stokes q and u of the standard star, finding 
q = 0.0093 ± 0.0014 and u = -0.0044 ± 0.0014. As this source has 
zero polarisation, these values can be seen as values for a compo- 
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Figure 2. This plot shows Stokes parameters q and u of objects in the com- 
bined rotl and rot2 datasets: bright field stars of the two datasets in green 
and blue (shifted to the (q, u) = (0, 0) point) and Swift J164449.3+573451 
in red. 



nent of the instrumental polarisation, and the Stokes parameters of 
the transient (and field stars) can be corrected with these. However, 
we first need to bring the observed values for the objects in the rotl 
and rot2 data to the same coordinate system as the standard star, 
through a simple rotation 

/ q' \ _ I cos 2(j> sin 20 \ I q 
\ u' J \ - sin 2(p cos 2<f> j\ u 

We note that this is not a complete calibration, as the standard and 
science data have a 90 degree rotation, i.e. the ordinary and ex- 
traordinary beams coming out of the Wollaston are swapped, so we 
can expect the effects of inefficiencies in the Wollaston (generally 
speaking of order P ~ 1%) to still be present after this calibration. 
However, the net polarisation of the transient can be well charac- 
terised through the foreground stars in the field of view. 

Three bright foreground stars (with q and u uncertainties 
smaller than 1%) are detected in both rotl and rot2 sets, and there- 
fore serve as the most useful secondary calibrators. There are sev- 
eral more sources detected in the frames, but we discard them based 
on proximity to the edges, noticeable asymmetry of the PSF or be- 
cause they are affected by bad pixels or strong vignetting. The av- 
erage of the field star q and ;/ values is somewhat offset from the 
(q, u) = (0, 0) point. This is the combined effect of polarisation in- 
duced by Galactic dust scattering (Galactic interstelar polarisation, 
GIP) and systematic errors (among others the non-ideal nature of 
the Wollaston prism, which can in an ideal situation be corrected 
for by using more instrument rotation angles), as we can assume 
that the average net intrinsic polarisation of field stars is zero. We 
compute the shift in q and u by taking the average of the field star 
values, and use the standard deviation as a conservative error. We 
then shift the sources in the q, u plane by this amount, and com- 
bine the rotl and rot2 values of the transient by taking the average. 
To the resulting error we add in quadrature the uncertainty in the 
average field star q, u value used for the shift. For the transient we 
find q = -0.035 ± 0.025 and u = 0.074 ± 0.021, i.e. Stokes U is 
non-zero with more than 3cr confidence. We plot the resulting q, u 
diagram in Figure [2] From this figure it can be seen that the tran- 
sient is significantly offset from the field stars, as expected of a truly 
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Figure 3. The top panel shows our WSRT radio lightcurves at 4.8 GHz 
and 1.4 GHz (red and cyan, respectively), in combination with EVLA 4.9 
GHz data from Zauderer et al. (201 1 ; blue circles). Note that in many cases 
the errorbars are smaller than the symbol size. The panel below that shows 
the 3cr limits on the linear polarisation at 4.8 GHz (red, left vertical scale) 
and 1.4 GHz (cyan, right vertical scale), with added the VLBA 8.4 GHz 
limits (open symbols, Levan et al. 201 1). The third panel shows the infrared 
lightcurves in J, H and K bands (red, green and blue squares, left vertical 
scale; photometry is not host subtracted), combining data reported in Levan 
et al. (2011) and Burrows et al. (2011) and additional data from Levan et 
al. in prep. The red star shows the WHT Ks band polarimetry (right vertical 
scale). The bottom panel shows the Swift XRT X-ray lightcurve, from the 
online Swift XRT lightcurve repository (Evans et al. 2007). All lightcurves 
are plotted with the BAT trigger time (12:57:45 UT on 28 March 2011; 
Burrows et al. 201 1) as zero time. 
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polarised source, but the uncertainties are substantial, as expected 
for a source this faint. 

Theoretical models are generally expressed in terms of the lin- 
ear and circular polarisation fractions P liu and P cil , and we therefore 
convert the Stokes parameter information to Py m using the equa- 
tions above, using the Stokes parameters q C0I ,u C0I : the Stokes pa- 
rameters obtained when we shift the average of the field stars to 
(q, u) = (0,0). Shifting the distribution by the field star average 
also eliminates the influence of polarisation induced by GIP, which 
in any case is expected to be low: with E(B - V)Gaiactic = 0.02 
(Schlegel et al. 1998) we expect this component of the polarisation 
to be < 0.5%. 

The uncertainty on the linear polarisation angle is a function of 
the intrinsic polarisation degree, cr g = a Plin l2P\ m , i.e. for low polar- 
isation values and faint fluxes the uncertainties in 6 are large. Errors 
on q and u are generally distributed as a Normal distribution and the 
Stokes parameters can have positive and negative values. In con- 
trast, Pjjj, is a positive definite quantity, and directly using the equa- 
tions above will therefore lead to an overestimated Py m and incor- 
rect confidence intervals, an effect generally referred to as the linear 
polarisation bias. The correction to the resulting and associated 
confidence ranges has been studied through both analytical and nu- 
merical (Monte Carlo) methods. Generally speaking, this correc- 
tion depends on o>, P, and the signal to noise (SNR) of /,. The 
Wardle & Kronberg (1974) prescription is often used in the litera- 
ture, in which the input P values are multiplied by - (a P /P) 2 
to find the bias-corrected polarisation. We follow Sparks & Axon 
(1999) in using a parameter q = P ■ SNR(/,) to trace the expected 
behaviour of the bias and cr F : when q > 2 the Wardle & Kron- 
berg correction is valid, and o> is as computed directly from the 
uncertainty of / ; . For the transient we find q = 2.3. 

From the relations above we find a formal measurement of the 
linear polarisation of the transient in the Ks band Pij n = 7.4 ± 3.5% 
and angle 9 = -32 ±12 degrees (corrected for polarisation bias). 

2.2 Radio 4.8 GHz polarimetry 

Radio lightcurves for Swift J164449. 3+573451 have been reported 
in Zauderer et al. (201 1) and Levan et al. (201 1), which can be re- 
produced by the emission from a relativistic jet. However, the only 
reported linear polarimetry of this source are from two epochs of 
VLBA observations at 8.4 GHz, giving 2<r upper limits of 4.5% 
and 4.7% (Levan et al. 2011); we plot these in Fig. [3] We report 
in this paper on radio observations performed with the Westerbork 
Synthesis Radio Telescope (WSRT) at 1.4 a nd 4.8 GHz. We used 
the Multi Frequency Front Ends dTarjl 199lh in combination with 
the IVC+DZB back end in continuum mode, with a bandwidth of 
8 x 20 MHz. Gain and phase calibrations were performed with the 
calibrator 3C 286 at both observing frequencies. The observations 
have been analysed using the Multichannel Image Reconstruction 
Image Analysis and Display (MIRIAD; ISault. Teuben & W right 
[19950 software package. 

We observed the source at six epochs, covering April 1 to May 
25, of which two observations were taken at 4.8 GHz and the other 
four at 1.4 and 4.8 GHz. The log and results of these measurements 
are shown in Table [T] The flux densities of the first three epochs 
have also been presented in Levan et al. (2011). Besides the flux 
densities, we have determined the linear polarisation by generating 
Stokes Q and U maps. The source was not detected in those maps 
at all epochs, and Table [T] gives the resulting 3cr upper limits on 
the linear polarisation P tm . Due to the low flux density at 1.4 GHz, 
the limits at this frequency are not well constrained, while upper 
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Flux 


Pv„, 
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(2011) 


(hours) 
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( c /r) 


March 31.904 -April 1.393 


12 


4.8 


990±21 


<9.7 


April 3.897 - 4.375 


5.3 


4.8 


1573±28 


<8.5 


April 9.880 - 10.359 


5.3 


4.8 


2185±30 


<6.0 


April 16.861 - 17.349 


12 


4.8 


2893±20 


<3.0 


May 6.807 - 7.285 


5.3 


4.8 


5877±28 


<2.1 


May 24.757 - 25.236 


5.3 


4.8 


6209±32 


<2.1 


April 3.920 - 4.395 


5.3 


1.4 


221 ±43 


<75 


April 9.903 - 10.379 


5.3 


1.4 


284±88 


<100 


May 6.830 - 7.305 


5.3 


1.4 


523±78 


<23 


May 24.780 - 25.256 


5.3 


1.4 


614±77 


<22 



Table 1. WSRT radio flux densities and 3cr upper limits on the linear polar- 
isation. 



limits at 4.8 GHz vary from ~ 10% on April 1 down to ~ 2% in the 
measurements in May 201 1. All limits and the resulting lightcurves 
are shown in Fig.fS] Both 4.8 GHz and 1.4 GHz show a clear, likely 
achromatic, lig htcurve break, w hich has been clearly seen at higher 
frequencies in Zaude rer! d201ll) and discussed in further detail in 
iMetzger. Giannios & Mimical < l201ll) . 



3 DISCUSSION 
3.1 Caveats 

There are a few caveats to point out. Firstly, there is considerable 
uncertainty on the fraction of the received K band flux that is pro- 
duced by the transient source, and what fraction is (stellar) emis- 
sion from the parent galaxy. The lightcurve monitoring and optical 
spectroscopy presented in Levan et al. (2011) shows that the tran- 
sient is very red, making the Ks band the most suitable optical/nIR 
wavelength range to maximise the transient vs. host galaxy light 
ratio. It also shows that the host level is likely significantly below 
the brightness of the source at the time of the polarimetry, which 
is further confirmed by late time infrared monitoring which shows 
that the source continued to fade significantly (at least 0.5 magni- 
tudes) at later times after the 80 day period shown here (Levan et 
al. in prep.; Fig. [3}- In addition, optical spectros copy shows brigh t 
emission lines and a fairly blue host continuum jLevanetal.l201ll) . 
As such, we believe that the majority of the light received in the Ks 
band is from the transient source. 

Secondly, the red colour of the transient likely implies consid- 
erable dust extinction within the host. Though the exact reddening 
A v depends on the assumed model used to fit the spectral energy 
distribution, it is l ikely to be in the range A v ~ 3 - 10 magnitudes 
(Bloo met alj201lb . Scattering of light onto dust grains induces po- 
larisation, whose orientation and magnitude as a function of wave- 
length depends on the geometry of the dust cloud with respect to the 
emitter and the dust grain size distribution. The high dust extinction 
values may correspond to a host galaxy dust induced linear polar- 
isation (HGIP) of a few percent. The HGIP should be vectorally 
added to the intrinsic polarisation (i.e. it adds in Q, U space, so it 
may increase the measured polarisation P lm or decrease it, depend- 
ing on 8). To derive a crude es t imate , we take the parametrisation 
from Mar tin. Clayton & Wolfll Jl999l) to describe p(/l)/p ma x, take 
Pmax/E(B - V) < 9%/mag and assume that all observed Ks band 
polarisation is caused by dust in the host, we find E(B - V) ~ 4, 
which appears too high compared to SED fits. Note that this also 
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implies a polarisation in optical bands a factor 3-4 higher than in the 
Ks band. However, it is clearly not possible to quantify the HGIP 
contribution: we have no information on the dust grain size distri- 
bution (as traced to some degree by R v ), the exact extinction value 
is very uncertain, and the applicability of Serkowski's parametrisa- 
tion of the relation between extinction and polarisation (Serkowski 
1973) is unclear. A firmer understanding of the extinction within 
the line of sight in the host galaxy is required before strong con- 
clusions can be drawn from the absolute value of the presented 
polarimetry. In particular, polarimetry at multiple broadband filters 
(or spectropolarimetry) can in principle differentiate dust scattering 
geometries (the received infrared light of Swift J 164449.3+573451 
may in fact be an echo, pure scattering) and dependence on grain 
size distribution (e.g. Zubko & Laor 2000). 

However, we note that the recently discovered source 
Swift J2058.4+051 6 shows remark able similarity to Swift 
J 164449.3+573451 dCenko et al]|201ll) . but has very low internal 
dust extinction, and shows a significant detection of linear polar- 
isation at very similar levels to this source at optical wavelengths 
(Levan et al. in prep.). Based on this similarity we proceed to con- 
sider the consequences of the bulk of the detected linear polarisa- 
tion in the Ks band being intrinsic to the emitting source. 

3.2 The nature of the radio and near-infrared emission and 
the expected polarisation 

Fits to the spectral energy distribution (SED) and time variable be- 
haviour of Swift J164449. 3+573451 performed at several differ- 
ent time slices and performed by different groups with indepen- 
dent data (Bloom et al. 2011; Zauderer et al. 2011; Burrows et 
al. 2011) broadly agree on a number of key points: the SED con- 
sists of two prominent bumps, peaking roughly at millimetre and 
X-ray/gamma-ray frequencies; the X-ray and radio emission are 
formed in different regions; the presence of a collimated, relativis- 
tic jet; the synchrotron nature of the long wavelength emission. An 
acceptable fit to the SED can be achieved with a model where the 
longer wavelength photons (i.e radio, (sub)mm) are produced by 
synchrotron emission, and the higher energy photons require the 
addition of synchrotron self-Compton (SSC) effects or further ex- 
ternal Compton (EC) processes (as in blazars). A somewhat poorer 
fit can be obtained assum i ng two independent s ynchrotron compo- 
nents ( iBloom et alj|201 j| ; iBurrows et al.ll20flh . requiring a larger 
extinction at optical wavelengths. The origin of the inverse Comp- 
ton seed photons i n the two component model is not clear (see e.g. 
IBloom et alj201 1|) , but disk photons seem a likely source. The ori- 
gin of the long wavelength synchrotron photons is more easily diag- 
nosed: the radio-microwave lightcurves and spectra are well repro- 
duced by the emission from the shock produced by th e interaction 
of a relativistic jet with the circumnuclear medium dBloom et al.l 



1201 ll ; IZauderej201 ll ; lMetzger. Giannios & Mimical201 ll) . The lat 
ter paper explores the emission from this shock in more detail, ex- 
plaining the achromatic break in the radio light curves (visible in 
Fig-HJ by the transition from a phase in which reverse shock emis- 
sion plays a role in the dynamics, to a phase where a Blandford- 
McKee self-similar evolution is approached. Our lightcurves show 
the first detection of this break at 1.4 GHz. 

These models, in which the long wavelength emission comes 
from an afterglow rather than in ternal synchrotron production 
within the jet (i.e. like a bla zar; Ivan Velzen. Kording & Falckd 
1201 ll ; iMiller & Gultekinlboi ll) . make it possible to use inferences 
from modelling of GRB blastwaves to understand the expected po- 
larisation properties. GRB afterglow polarimetry at optical wave- 



lengths generally shows low but variable linear polarisation levels 
(below a fe w percent), and no indic ations of a significant circular 
polarisation dGranot & Tavlodl2005l ; Wiersema et al. in prep.). The 
observed variability in linear polarisation is often seen to track vari- 
ations in afterglow brightness on top of the g eneral power law de - 
cay (most dramatically seen in GRB 030329, Greine r et al. I l2003h . 
This variability is often explained in a scenario where the received 
emission from a surface of equal arrival time (GRB blastwaves are 
highly relativistic) consists of a large number of patches, each with 
a fairly coherent field and high polarisation. As the received emis- 
sion is the sum over many such patches, th e observed linear po- 
larisation is low (Gruzinov & Waxman|[l999h . In the case of Swift 
J164449. 3+573451, an observed polarisation of ~ 7% would imply 
~ 100 such patches (assuming a ~ 70% polarisation of synchrotron 
emission with a perfectly ordered field). Bumps in the GRB after- 
glow lightcurve are formed when distinct bright patches are mo- 
mentarily dominating the received emission, thereby dominating 
the received polarisation too. Radio polarimetry of GRB afterglows 
is highly challenging because of the low fluxes, but predictions in- 
dicate detectable levels of polarisation dToma. Ioka & Nakamural 
2008), with similar polarisation percentages to the optical polar- 
isation (but note that nearly all GRB afterglows are too faint at 
radio wavelengths to perform polarimetry at percent level). The 
non detection of radio polarisation at any epoch in the later time 
WSRT data is in agreement with the expectation from models ex- 
plaining the source SED and lightcurves with an afterglow-like 
phenomenon. The upper limits from WSRT can be used to probe 
the physics of the jets in detail o nce more of the microphysical pa- 
rameters can be derived (see e.g. lToma. Ioka & Nakamural2008l for 
details). 

lBloometal.N201ll) explicitly indicate that the expected lin- 
ear polarisation of this source should be relatively low, i.e. like a 
GRB afterglow. As explained above, we confirm this prediction at 
radio w avelengths. The modelling of Metzg er. Giannios & Mimica 
a201 ll) of a larger dataset using this scenario invokes the pres- 
ence of a reverse shock: the early rise in the radio light curve is 
occurring while the reverse shock is still passing through the ini- 
tial shell of ejecta. However, similar to GRB afterglows, the emis- 
sion is dominated the foward shock (i.e. the reverse shock only af- 
fects the forward shock emission through the dynamics of the jet 
expansion). Reverse shock emission is occasion ally detectable in 
GRB afterglows (bu t appears to be fairl y rare; iKlotz et alJ|2009l; 
iMelandri et al]|2008t iRvkoff et alj|2009l) . and may be highly po- 
larised, in contrast to forward shock emission, for example in a 
scenario where there is both a large scale, weak, ordered mag- 
netic field in the ambient medium and a tangled, random, fiel d 
generated by postshock turbulence (e.g. iGranot & K6niglll2003l) . 
In this case the polarisation P varies as a result of changes in 
the ratio of the ordered- to-random mean-squared field amplitudes 
dGranot & K6niglll2003l) . This model was put forward by the au- 
thors to explain the absence of large variation in P and a rota- 
tion of 9 at the time of the so-called jet-break, expected if GRB 
afterglow emission is synchrotron emission domin ated by a tan- 
gled magnetic field viewed somewhat off-axis (e.g. lLazzatill2006l 
and references therein). We detect no polarisation at radio wave- 
lengths during the steep rise phase in the Swift J164449. 3+573451 
lightcurves with fairly deep limits (Fig. [3 ] >, whic h agrees with the 
findings of iMetzger. Giannios & Mimical d201 ll) that the received 
emission is dominated by forward shock emission. The non de- 
tection of any radio polarisation at any epoch also excludes large 
amplitude variability of P, and is consistent with the absence of a 
jet-break like phenomenon (or other geometric, symmetry break- 
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ing effect s) under the assumption of a pure tan gled magnetic field. 
Note that lMetzger. Giannios & Mimical fcoill) infer that 8 ja < l/T 
throughout: the entire jet is observed, similar to GRBs observed 
long after the jet break. Further detailed modelling of the jet dy- 
namics will be required to determine the influence of this configu- 
ration on the predicted polarisation signal. 

Interestingly, the infrared wavelengths find themselves in the 
region where there are contributions both from the synchrotron 
component seen in radio wavelengths, and from the second bump 
in the SED, where inverse Compton effects play an important role. 
The detection of a low level of polarisation, described above, is not 
inconsistent with the expectations from an afterglow origin of the 
long wavelength emisison, and is comparable with GRB afterglows 
(see Wiersema et al. in prep, for Ks band polarimetry of the after- 
glow of GRB 091018), though 7% polarisation is high for a GRB 
afterglow. This value is in l ine with the non-detection of ra dio po- 
larisation in this scenario dToma. Ioka & Nakamural l2008). How- 
ever, there is a possibility for a contribution of non-synchrotron 
emission processes to the measured polarisation. The expected po- 
larisatio n properties of invers e Compton emission have been eval- 
uated bv lKrawczvnskil ilOl lh in the framework of blazar emission, 
who demonstrates that SSC emission from unpolarised seed pho- 
tons produces vanishing linear polarisation for electron Lorentz 
factors y > 10. While the origin of the seed photons for this 
source is not well established, it see ms likely that they originate 
from somewhere in the accretion disk l lBloom etalj201lh or are the 
same photons that form the radio-optical spectrum (see lAliu et al.l 
l201ll for a discussion), so in both cases are effectively unpol arised. 
The analytical and numerical results o f lKrawczvnsklbOllh indi- 
cate that it is unlikely that SSC contributes to the observed Ks band 
polarisation. 

We can now contrast the polarisation properties of Swift 
J 164449.3+573451 with other transient relativistic jet sources, in 
particular X-ray binaries. In these sources the detected synchrotron 
emission is internal to the jet, not from an afterglow like phe- 
nomenon. As such, their polarisation properties are different in 
some aspects from GRB afterglows. Radio polarimetry of low mass 
X-ray binary jets shows the polarisation and its position angle to 
be a function of wavelength and spectrum. For optically thick jets 
(flat spectrum, self-absorbed) the synchrotron polarisation is low 
(a few percent ), whereas fo r optically thin jets it is much higher 
(10 - 30%; see lFendeJl2006l for a review). In the optical and near- 
infrared wavelengths there may be a contribution from Rayleigh 
scattering in the vicinity of the accretion disk (the massive star 
plays a role as well in high mass X-ray binaries), and when jets 
are seen these wavelength s diagnose the magne tic field closer to 
the base of the jet (e.g. iRussell & Fen der 2008). Sometimes this 
is optically thin and the pol arisation is then seen to be a few per- 
cent and highly variable (e.g. lRussell & Fender 2008). Of particular 
interest for comparison with Swift J 164449.3+573451 is the low 
mass X-ray binary (black hole candidate) XTE J1550-564: during 
a faint X-ray ou tburst in the low-hard state, Ks band polarimetry 
was obtained by Dubu s & Chatvl ( 12006] ) . who find an intrinsic (i.e. 
after correcting for the Galactic dust induced polarisation through 
field stars) linear polarisation in the Ks band of 0.9% - 2% (95 % 
confidence interval). The authors attribute this polarisation to op- 
tically thick synchrotron emission from a compact jet of the bi- 
nary. Their measurements of low level polarisation in this case, 
far from the self-absorption frequency, is comparable to the polar- 
isation of Swift J164449.3+573451. When the origin of the near- 
infrared emission of Swift J 164449.3+573451 is better understood 
(through long timescale lightcurves and late time SEDs), the com- 



parison with XTE J 1550-564 and X-ray binaries in general can be 
further quantified. 

Finally, the comparison of Swift J164449.3+573451 with 
blazars has been made seve r al tim es in the literature. As 
also remarked in iLevan et"aT, the polarisation of this 

source is markedly lower than sometimes observed in blazars 
dAller. Aller & Hughes 2003) but it is generally consistent with 
what is seen when they are in a relatively low flux state: the most 
strongly optical / near-infrared polarised emission appears to corre- 
spond to the peak of blazar flares. Our LIRIS observations of Swift 
J 164449.3+573451 occurred at relatively late times, and the evo- 
lution of the event was much more rapid than typical blazar flares. 
Multi-epoch polarimetry of events like Swift J164449. 3+573451 
when they are closer to the peak brightness may therefore shed 
more light on the link of this source with blazar polarisation vari- 
ability. 

As a last point we wish to point out the suitability of sources 
like Swift J164449. 3+573451 for X-ray polarimetry missions, for 
example GEMS, which would have been able to perform percent- 
level polarimetry in a reasonable exposure time (~ 100 ks) at a few 
days after trigger (see Fig. [3]for the X-ray lightcurve). The X-ray 
and radio brightness of Swift J 164449.3+573451 and its transience 
(i.e. tracking the synchrotron break frequencies) allows for new in- 
sights into the physics of jet launching, particle acceleration and 
photon emission processes: X-ray polarimetry can distinguish be- 
tween external Compton and SSC effects, while ground-based (ra- 
dio - optical) studies can determine the micro- and macrophysics of 
the jet-medium interaction. 



4 CONCLUSIONS 

We present in this paper a measurement of the Ks band linear po- 
larisation of Swift J 164449.3+573451, giving Stokes parameters 
Q/I = -0.035±0.025 and U/I = 0.074±0.021, or P hn = 7.4±3.5%. 
In addition we present an extensive set of upper limits on the lin- 
ear polarisation at two radio frequencies, 4.8 and 1.4 GHz. These 
measurements confirm the results from SED modelling that the 
long wavelength emission originates in the interaction of a rela- 
tivistic jet with the medium around the black hole in the centre 
of the host galaxy. The SSC component inferred from SED mod- 
elling is not likely to contribute to the observed polarisation in the 
Ks band. It seems likely that a polarisation component caused by 
HGIP is present in the Ks band polarimetry. It would require a bet- 
ter knowledge of the shape of the extinction curve and the amount 
of extinction A K to correct for this. However, the recent detection 
of Swift 2058.4+05 16 shows that relatively unobscured versions 
of Swift J164449.3+573451 do exist, where HGIP is negligible. 
In those less obscured cases multi-colour polarimetry can be ob- 
tained, which allows for HGIP correction in a direct way. Our ra- 
dio polarimetry restricts the possible configuration of the magnetic 
field and its coherence through limits on polarisation during both 
the steep rise phase in the radio lightcurves and deep limits at late 
times. Lastly, our study demonstrates the possibility to probe these 
very faint transients using 4m class telescopes (with highly versa- 
tile instruments like LIRIS). 



ACKNOWLEDGMENTS 

We thank the ING staff for their support of the LIRIS polarime- 
try effort discussed here, in particular M. Hrudkova. We thank the 



Polarimetry ofGRB 110328 /Swift J164449.3+573451 7 



anonymous referee for useful suggestions. KW acknowledges sup- 
port from STFC. RLCS is supported by a Royal Society Fellow- 
ship. The William Herschel Telescope is operated on the island of 
La Palma by the Isaac Newton Group in the Spanish Observatorio 
del Roque de los Muchachos of the Instituto de Astrofsica de Ca- 
narias. The WSRT is operated by ASTRON (Netherlands Institute 
for Radio Astronomy) with support from the Netherlands founda- 
tion for Scientific Research. This work made use of data supplied 
by the UK Swift Science Data Centre at the University of Leicester. 



Wardle, J. F. C. & Kronberg, P. P. 1974, ApJ, 194, 249 
Woosley, S. E. & Heger, A. 2011, ApJ Submitted, 
larXiv: 11 10.38421 

Zauderer, B. A. et al. 2011, Nature, 476, 425 
Zubko, V. G. & Laor, A. 2000, ApJS, 128, 245 



REFERENCES 

Aliu, E. et al. 201 1, ApJ, 738, 30 
M. F. Aller, H. D. Aller, P. A. Hughes 2003, ApJ 586, 33 
Bloom, J. S. et al. 201 1, Science 333, 203 
Bower, G., Bloom, J., & Cenko, B. 2011, ATel, 3278, 1 
Burrows, D. N. et al. 201 1, Nature, 476, 421 
Cannizzo, J. K., Troja, E. & Lodato, G. 2011, ApJ, 742, 32 
Cenko, S. B. et al. 2011, ApJ submitted. larXiv: 1 107.53071 
Chandrasekhar, S. 1960, Radiative transfer, New York: Dover, 
1960 

Cummings, J. R. et al. 2011, GCN 11823 
Dubus, G. & Chaty, S. 2006, A&A, 458, 591 
Evans, P. A., et al. 2007, A&A, 469, 379 

Fender R. P., 2006, in Lewin W. H. G., van der Klis M., eds, Com- 
pact Stellar X-Ray Sources. Cambridge Univ. Press, Cambridge, 
p. 381 

Granot, J. & Konigl, A. 2003, ApJ, 594, L83 
Granot, J. & Taylor, G. B. 2005, ApJ 625, 263 
Greiner, J. et al. 2003, Nature, 426, 157 
Gruzinov, A., & Waxman, E. 1999, ApJ, 511, 852 
Klotz, A., Boer, M., Atteia, J. L., & Gendre, B. 2009, AJ, 137, 
4100 

Krawczynski, H. 2011, ApJ in press. larXiv: 1109.21861 
Krolik, J. H. & Piran, T. 201 1, ApJ, 743, 134 
Lazzati, D. 2006, NJPh, 8, 131 
Levan, A. J. et al. 2011, Science, 333, 199 

Manchado, A., Barreto, M., Acosta-Pulido, J.A., et al. 2004, in 

Proc. of the SPIE, 5492, 1094 
Martin, P. G., Clayton, G. C, Wolff, M. J. 1999, ApJ, 510, 905 
Melandri, A., et al. 2008, ApJ, 686, 1209 

Metzger, B. D., Giannios, D. & Mimica, P. 2011, MNRAS sub- 
mitted, larXiv: 11 10. lTTTI 
Miller, J. M. & Gultekin, K. 201 1, ApJ, 738, 13 
Oliva, E. 1997 A&AS, 123, 589 

Quataert, E. & Kasen, D. 2011, MNRAS in press, 
larXiv: 1105.32091 

Russell, D. M. & Fender, R. P. 2008, MNRAS, 387, 713 
Rykoff, E. S. et al. 2009, ApJ, 702, 489 

Sault, R.J., Teuben, P.J., & Wright, M.C.H. 1995, Astronomical 

Data Analysis Software and Systems IV, 77, 433 
Schlegel, D. J., Finkbeiner, D. P. & Davis, M. 1998, ApJ, 500, 525 
Serkowski, K. 1973, in Interstellar Dust and Related Topics, ed. 
J. M. Greenberg, & H. C. van de Hulst (Dordrecht: Reidel), IAU 
Symp., 52, 145 

Shao, L„ Zhang, F.-W., Fan, Y.-Z., Wei, D.-M. 201 1, ApJ, 734, 33 

Sparks, W. B. & Axon, D. J. 1999, PASP 111, 1298 

Tan, G.H. 1991, in ASP Conf. Ser. 19: IAU Colloq. 131: Radio 

Interferometry. Theory, Techniques, and Applications, 42-A6 
Toma, K., Ioka, K. & Nakamura, T. 2008, ApJL, 673, 123 
van Velzen, S., Kording, E. & Falcke, H. 2011, MNRAS, 417, 51 



